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The effects of the axial anomaly are suppressed at high temperatures due to screening effects 
in the quark-gluon plasma. If the suppression is nearly complete close to the chiral transition 
temperature, this can have consequences for the nature of the phase transition. The use of a chiral 
action such as Domain Wall Fermions allows us to gain a deeper insight into the issue. Our lattice 
sizes were 16^ x 8 x L,, with Lj = 32 or 48, and our pion mass was approximately 200 MeV. We 
found that t/ ( 1 )a stayed broken above the chiral transition. However the breaking was found to be 
due to topologically nontrivial configurations which raises the question as to whether it persists in 
the thermodynamic limit. We also present results for the eigenvalue density of the Dirac operator 
It is seen that although the density decreases dramatically across the chiral transition temperature, 
U{1)a still remains broken at our current volume and quark mass due to the presence of zero 
modes. 
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1. Introduction 

The Lagrangian of Quantum Chromodynamics (QCD) with Nf massless flavors of quarks 
is invariant under a global SU {Nf)L (X" SU (A^/)r 0U{l)v ® U{1)a symmetry. In the vacuum, the 
SU {Nf)L(!$SU {Nf)R chiral symmetry is spontaneously broken to a SU {Nf)v subgroup, correspond- 
ing to flavor symmetry. This spontaneous breaking of chiral symmetry gives rise to a nonvanishing 
expectation value (yy) of the chiral condensate. 

The axial U{1)a symmetry of the QCD Lagrangian on the other hand is broken by the axial 
anomaly. The inclusion of quantum fluctuations leads, at the perturbative level itself, to non- 
conservation of the axial current; this is the famous Adler-Bell-Jackiw anomaly [[Tj, 

{d,n = -^{^"^''V;s)- (1-1) 

In QCD, the anomaly implies global non-conservation of axial charge. Naively, integrating 
Eq. ( |1.1[ ) over all spacetime should give zero since the left-hand side is a total divergence. However 
there exist special gauge field configurations in QCD for which the integral of the right-hand side 
is not zero. These are the configurations with nontrivial topology [||]. All such configurations must 
be included in the path-integral. Anomalous contributions arise for any observable for which the 
contribution from such configurations is unsuppressed. 

1.1 Effective U(1)a Restoration 

A common example of a phase transition in several finite-temperature field theories is the 
restoration of a spontaneously broken global symmetry. This is the case with chiral symmetry in 
QCD as well. For Nf = 2 in the massless limit, the phase transition is expected to be second-order 
and belonging to the 0(4) universality class. When the quarks are massive, this transition becomes 
a crossover. 

By contrast axial symmetry is broken at the perturbative level itself. There is thus no question 
of its complete restoration at any temperature. However as we have already seen, anomaly-related 
effects arise from the existence of topologically nontrivial configurations. The action for these 
configurations is proportional to cc^^. Such actions are therefore Boltzmann-suppressed due to 
the screening of the coupling constant at high temperatures [Q]. Although there is always some 
amount of U{1)a breaking below T = oo, it is conceivable that this suppression is nearly complete 
by some temperature that is not too high. We may then speak of an effective restoration of the axial 
symmetry. 

If this temperature is close to the chiral phase transition temperature T^, then the effective 
restoration of ?7(1)a can have interesting phenomenological consequences. The standai^d picture 
of a second-order phase transition is based on the assumption that U{1)a breaking is substantial 
near r^.' If this is not the case, then the phase transition may even be first order [||]. Understanding 
the contribution of ?7(1)a is thus essential to mapping the phase diagram of QCD. 

'Note that the chiral condensate (yv)' which signals chiral symmetry breaking, also breaks U{1)a. Consequently 
there is no question of U(1)a being restored before SU {2)l x SU (2)r is. 
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2. Domain- Wall Fermions and DSDR 

Chiral symmetry restoration and effective U{\)a restoration at high temperatures are both 
non-perturbative phenomena whose rehable study demands the use of nonperturbative techniques. 
Currently, lattice QCD is certainly the most viable and reliable such technique. Extensive lattice 
QCD studies of chiral symmetry restoration have already been carried out (for a review and sum- 
mary see [§, ^). The question of ?7(1)a restoration too has been investigated before ^ [TTl ]. 



However such studies have almost always been carried out with staggered fermions. For these 
fermions the issues of chiral symmetry, the anomaly and the relation between the anomaly and the 
index theorem are very subtle [|T2|, [13], Q. Hence further studies using different fermion discretiza- 
tion schemes are certainly welcome. 

Domain Wall Fermions are a fermion discretization scheme that preserves the full SU {Nf )i x 
SU{Nf)R chiral symmetry of continuum QCD and also reproduces the correct anomaly even at 
nonzero values of the lattice spacing |jT5|]. The domain wall formulation is one of five-dimensional 
fermions whose low-energy spectrum is four-dimensional and also, when the fifth dimension is in- 
finite in extent, exactly chiral. The gauge fields remain four-dimensional and couple to the fermions 
in the usual way. For finite fifth dimension, the residual chiral symmetry breaking manifests itself 
at low energies as an additive shift mpes of the bare quark mass 



The QCD phase transition has been studied before with domain wall fermions [ ]17| , |18[ ]. A 
challenge encountered in the most recent study was the rapid variation of wires as one moved to- 
ward stronger coupling which made it difficult to keep the pion mass fixed throughout the tem- 



perature range studied [18]. The use of improved gauge actions such as the Iwasaki action results 
in a smaller value of mres overall but cannot arrest the rapid growth of wires as the temperature is 
decreased. 

In an ongoing study of QCD thermodynamics using domain wall fermions by the HotQCD 



collaboration [19], mres was sought to be kept to a minimum through the use of the "Dislocation 
Suppressing Determinant Ratio (DSDR)." The usual Iwasaki gauge action was augmented with 
a ratio of Wilson determinants which suppressed the zero modes (dislocations) that contributed 
to mres- To maintain adequate topological tunneling, the Wilson-Dirac mass was set equal to the 
domain- wall height — Mq plus a small chirally twisted mass ie/j viz. [ pO| , ^ 



det 



det 



Dlf ( -Mo + iehY5 )Dw{-Mo + istYs ] 



(2.1) 



With this action we generated a few thousand configurations each at seven temperatures be- 
tween 140 MeV and 200 MeV. Our lattice sizes were 16^ x 8 x L, with L, = 32 for T ^ 160 MeV 
and Ls = 48 at lower temperatures. Before generating these configurations, we generated several 
zero-temperature ensembles at several values of the coupling j8 both to set the scale and to deter- 
mine the residual mass mres- The input light and strange quark masses were then chosen so as to 
keep the kaon physical and the pion mass fixed at 200 MeV for all j8; this defined our Line of 
Constant Physics. 

At each temperature, we measured (^^f) , ("V?75 ^) and the corresponding disconnected suscep- 
tibiUties. We also measured the flavored scalar (5), pseudoscalar (tt), vector (p) and axial vector 
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(fli) correlators to be discussed below. Separately, we also measured the topological charge 2top 
for each configuration through the use of cooling and smeared gauge field operators [23]. From 
this we calculated (2top)> (|2top|) and the topological susceptibility ;^top- Finally, we also measured 
the lowest hundred eigenvalues of the five-dimensional Dirac operator on each configuration in an 
effort to determine the eigenvalue density distribution p(A) (Section |5]). The physics behind the 
DSDR action, its performance and the results for the chiral phase transition have been presented by 
M. Cheng at this conference [24]. A complete description of our ensembles, scale determination 
and measurements is also forthcoming [p^]. 
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Figure 1: The one-flavor disconnected chiral susceptibility for the light quark. The transition region is 
broad, as might be expected of a crossover, with a peak near T « 160 MeV. 



Fig. [T] plots the disconnected chiral susceptibility as a function of the temperature. The sus- 
ceptibility peaks at 160 MeV; accordingly we take that to be the approximate value of the chiral 
phase transition temperature Tc. Since the phase transition is expected to merely be a crossover for 
ruji > 0, this value only serves as a reference when discussing the possibility of ?7(1)a restoration. 



3. Symmetries, Correlators and Susceptibilities 



The influence of a symmetry is seen on the appropriate correlators of Dirac bilinears. In the 
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Figure 2: The symmetries SU{Nf)L x SU{Nf )R and U{1)a relate mesons in different spin-flavor channels. 
The above diagram summarizes these relations for Nf = 2. 
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scalar-pseudoscalar sector, for Nf = 2, we have two iso-triplet correlators: 

Cs{x) = (ud{x) du{0)^, 

as well as two iso-singlet correlators: 



Caix) = [uu{x) +dd{x)) (mm(0) +dd{0)) 
Crj' W = ( (iM75M(x) +it/75(i(x)) (iM75M(0) +ic^75^/(0)) ^. 



(3.1a) 
(3.1b) 



(3.2a) 
(3.2b) 



The 5 and 71 correlators receive contributions only from diagrams with connected quark lines and 
are thus easier to measure. The a and the rj' on the other hand receive contributions from diagrams 
with connected as well as disconnected quark lines. The connected parts of these correlators are just 
the 5 and the n respectively. The full correlator however is obtained only after this part is canceled 
by a similar contribution from the disconnected piece. This makes the a and rj' correlators much 
harder to measure. 

These correlators transform into one another under a chiral or an axial rotation, as summarized 
in Fig. 1^ This implies for e.g. that the n and the a (5) correlators become identical when chiral 
(axial) symmetry is restored. 



In addition to the correlators in eq. (|3JJ) we also calculated the connected vector and axial 
vector correlators viz. 



Cp{x) = {uY^d{x) dY^u{0)), Q,(;c) = {iuysy^dix) idY57^lu{0)) 



(3.3) 



An axial rotation has no effect on the vector and axial vector correlators. The two are in 
fact related through chiral transformations and they become degenerate when chiral symmetry is 
restored. We plot these correlators at T = 150 MeV and 160 MeV in Fig. ||. Since the transition 
is a crossover the two only become exactly identical at very high temperatures, but they are nearly 
degenerate by T = 160 MeV. 
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Figure 3: The vector (p) and axial vector (ai) correlators for T = 150 MeV and 160 MeV respectively. 

By integrating these correlators over the four-volume, we obtain the corresponding suscepti- 
bilities Xn, Xa, etc. Just as for the correlators, one has connected and disconnected susceptibilities 
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depending on the type of correlator being integrated. Furthermore, the disconnected parts of the a 
and T]' susceptibilities are equal to the disconnected susceptibilities ;^disc and XsAi^c viz. 



' iW"^) f = Xdi^c and X77',disc = ( {WYswf ) = Zs.disc (3.4) 



The appropriate symmetry restoration gives rise to equalities among the different susceptibilities: 



Xn = Xs + X&^c and 
Xn = Xs and 



X5=Xn- Z5,disc. [SU{2)l X SU{1)r] 
XS + Zdisc =Xk- Z5,disc • [U{1)a]- 



(3.5a) 
(3.5b) 



The difference Xn — Xs must go to zero as ?7(1)^ breaking is suppressed. Eq. ( |3.5a| ) tells us that this 
difference equals ;^disc once chiral symmetry is restored. Moreover, we see that chiral symmetry 
restoration implies that XHfic = XsAisc whereas axial symmetry restoration implies the opposite, 
namely ;^disc = —Xs.disc- Either way, when both chiral and axial symmetry are restored, one has 
Zdisc = = X54isc- In other words, U (1)a restoration is signaled by a vanishing disconnected chiral 
susceptibility. 
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Figure 4: The susceptibilities Jdisc. Zs.disc and Xn ^ Xs for each of the temperatures. All are very nearly 
equal from T ^ 110 MeV onward. None of these susceptibilities vanishes for all the temperatures shown 
here. The red and blue points have been horizontally displaced by ±1 MeV for clarity. 

Fig. ^ plots these susceptibilities for each of the temperatures that we studied. Although the 
equalities derived in Eqs. (3^) are strictly valid only in the chiral limit, we see that ;(^disc> Z5,disc and 
Xtt — Xs are almost equal to each other from about 170 MeV onwards. Furthermore, none of these 
susceptibilities is equal to zero even atT = 200 MeV, the highest temperature that we studied. If 
we take Tc ~ 160 MeV, this would seem to suggest that U{1)a remains broken even at T ?a 1.257^. 



4. The Correlation with Topology 

Let us take a closer look at the source of U{1)a violation. If we write the n and 5 correlators 



(Eqs. (|3.lD) in terms of their left- and right-handed components, we get 



Cs/A^) = {uLdR{x)dRUL{0) + URdL{x)dLUR{Q)) 
± {uLdR{x)dLUR{0) + URdL{x)dRUL{0)) . 



(4.1) 



6 



Looking for U{1)a Restoration 



Prasad Hegde 



it + 6: 140-200 MeV 



t * f « ft 
I I <!> ^ 



10 12 14 



1- 5: 140-200 MeV 

' ^ in m -a * * 



10 12 14 



Figure 5: (Left) The sum of the 5 and n correlators. The temperature increases from 140 to 200 MeV as 
one moves downward along the y-axis. (Right) The difference of the two correlators. The temperatures are 
identified by the same symbols as in the plot on the left. 



Here the left- and right-handed parts are defined as 
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In terms of these, our scalar and pseudoscalar correlators are 

u{x)d{x) =ui{x)dR{x) + UR{x)di{x) and u{x)yid{x) = ui{x)dR{x) — UR{x)di{x) 
AU{\)a transformation is given by 



->e ^^ul{x), 

ur{x) 



t+'%R{x), 



ur{x) — )■ UR{x)e, 
ui{x) ^ul{x)&^'^ , 



(4.2) 



(4.3) 



(4.4) 



(4.5) 



and similarly for d{x). Looking back at Eq. ( |4. l[ ), we see that U{1)a violation comes entirely 
from the terms on the second Une, which occur with opposite signs for the correlators. By contrast 
the terms on the first line, which occur with the same sign for both correlators, are invariant with 
respect to ?7(1)a transformations. From this it is clear that the U{\)a violating and respecting parts 
may be isolated by looking at the 71 — 5 and the 71 + 5 correlators respectively. 

We plot the sum and the difference of these two correlators in Fig. |[ We see that the difference 
is of the same order of magnitude as the sum at the farthest separations (x ^ Na/1) for all the 
temperatures shown here. This reaffirms our earlier observation (Fig. |^ that U{\)a remains broken 
even at the highest temperatures that we studied. 
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Figure 6: Time histories of the integrated correlators (red lines) and the topological charge (blue lines) for 
T = 170-200 MeV. The topological charge histories have been displaced downward by 2-3 units for clarity. 



It is important to determine whether the U{1)a breaking that we observe is due to the presence 
of topologically nontrivial configurations or occurs merely because of the quark mass. The connec- 
tion to topology is confirmed when we look at the time histories of these correlators and compare 
these with the time histories of Qtop- Fig- 1 shows these time histories for four temperatures viz. 
170, 180, 190 and 200 MeV. To remove the dependence on a particular separation x, we plot the 
time histories of the integrated correlators viz. {7l{x) — 5{x)) = n — 5. In all but a few cases, 
the spikes in tt — 5 are found to line up with the jumps to nonzero values of the topology. In other 
words, U{1)a is broken not "on average" but rather by specific configurations whose frequency 
decreases as the temperature is increased. 



5. The Spectrum of the Dirac Operator 

The connection to topology is intriguing, but it also raises questions about the eventual fate of 
U{1)a breaking. Nontrivial topologies are distinguished by the fact that the Dirac operator always 
has a zero eigenvalue in their presence.^ The contribution of these modes however vanishes when 
the four-volume is sent to infinity i.e. in the thermodynamic limit. 

^The Atiyah-Singer theorem constrains the difference between the number of left- and right-handed zero modes viz. 
A'-i- — A'_ = 2top- Out of the total A'-i- +A'_ zero modes however, only 2top are stable with respect to small deformations. 
A stronger statement therefore is that in the presence of a configuration with winding number (gtop, the Dirac operator 
has 2top exact or robust zero modes. 
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To see this, let us express (v^V^) and X7t~X8 i^i terms of the eigenvalues A of the Dirac operator 

viz. 

(W)= rdAp(A,m)^^ + i^, (5.1a) 
Jo WJ + ^ mV 

n ^ 2(|etop|> 

The usual weighted average over gauge fields is re-expressed as an average over eigenvalues dis- 
tributed according to the spectral density p{X,m). The second term on each RHS represents the 
contributions of the exact zero modes [^. The first term represents the contributions coming from 
the rest of the spectrum. 

For large volumes, the topological charge is expected to obey a Gaussian distribution with a 
width proportional to the volume viz. |26] 

Eq. (|^ implies that (|2top|) °^ VV, hence the second terms in Eq. ( ^ ) vanish as V — )- oo. 

?7(l)A-breaking then must come from the rest of the spectrum i.e. from the first terms of 
Eqs. < \5.l\ ). In the chiral limit the dominant contribution to these integrals comes from the eigen- 
values within a small distance of the origin. This is similar to what happens when chiral symmetry 
is broken: Eigenvalues A ~ ^(l/V) i.e. the near-zero modes, build up near the origin and it is 
these, rather than the exact zero modes, that break chiral symmetry [p7|]. This is reflected in the 
Casher-Banks relation (yv) = ^P(OiO) for e.g. @]. 

The spectral density p{X,m), for small A, can be determined by looking at the distribution of 
the lowest eigenvalues of the Dirac operator with respect to the gauge configurations. For domain 
wall fermions, the correct Dirac operator is the four-dimensional one whose exact form unfortu- 
nately is unknown. However, since it is realized in the low-energy limit of the five-dimensional 
theory, its low-lying spectrum will be the same as that of the full five-dimensional theory upto 
an overall renormalization factor which may also be determined non-perturbatively. We show the 
resulting histograms for p(A) in Fig. ^. These results, as well as details of the renormalization 
procedure, have all been presented by Z. Lin at this conference [29] and will also be described in a 
forthcoming publication [p^. 

From Fig. ^ we see that the while the eigenvalue density at the origin shrinks dramatically in 
going across the chiral phase transition from 150 to 170 MeV, small eigenvalues still occur with 
reasonable frequency upto T = 180 MeV. On the other hand, for T = 190 and 200 MeV most of the 
eigenvalues occur away from the origin. Nevertheless at both temperatures there is also a second 
set of eigenvalues that occurs close to the origin. It is these eigenvalues that are responsible for 
U{\)a breaking. 

Fig. |8] shows the correlation between the magnitude of ?7(1)a breaking (i.e. the value of the 
integrated correlator n — 5) and the value of the smallest eigenvalue for the Dirac operator on a 
configuration-by-configuration basis. It is clear from the plot that the closer Amin is to zero, the 
larger the value of ;r — 5 for that configuration. Furthermore, Amin is well-separated from zero for 
the configurations with 2top = 0. Accordingly for such configurations tt — 5 is quite small. 
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Figure 7: (Left to right, top to bottom) The renormalized eigenvalue spectrum for T ~ 160 - 200 MeV. 
These figures have been taken from pi]. The histograms are w.rt. A, which in the continuum is related 



to the conventional Dirac eigenvalue X by A = ^J + (ini +OTres)^- The leftmost line in each plot marks 
the location of A = w; +OTres- In the continuum, the eigenvalue density at that point yields the value of 

P(0,OT, +mi.es)- 



The dominant contribution to Xn — Xs comes from the configurations with 2top 7^ 0. As ex- 
pected, Amin is quite close to zero for these configurations. These eigenvalues however shall vanish 
in the thermodynamic limit, and the magnitude of ?7(1)a breaking at T = 200 MeV and 190 MeV 
is likely to be much smaller than what we currently see. 

6. Summary and Discussion 

The question as to whether U{1)a is restored at very high temperatures has a long history. As 
already mentioned, such a restoration cannot be complete below T = oo. Nevertheless its breaking 
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Figure 8: Scatter plot showing the correlation between n — 5 (y-axis) and the value of the smallest eigenvalue 
Ajnin (x-axis) at r = 200 MeV on a configuration-by-configuration basis. A^m is related to the smallest 

Dirac eigenvalue by Amin — \J X^^^ + {mi + m^s)^- The dotted vertical lines bracket the unphysical region 
Amin ^ '«/ +'«res- Open circles denote configurations with Qtop = while bursts denote configurations with 
Stop 7^ 0. Configurations with large values of ;r — 5 are also those with 2top 7^ 0, indicating that U{1)a 
breaking at this temperature comes from the exact zero modes. Note that the y-axis is logarithmic. 

is very much suppressed at high temperatures and, depending on the magnitude of suppression, it 
could be effectively restored above some temperature. It is natural to ask whether this temperature 
is close to the familiar chiral transition temperature. In this work, we investigated this question in 
the context of the 2 + 1 -flavor theory, on the lattice, by working with a chiral action and examining 
the behavior of the scalar and pseudoscalar iso-triplet correlators. 

We found that U{1)a remained broken even after the usual SU {2)i x SU (2)^ chiral symmetry 
had been restored. Moreover this was due to the presence of configurations with 2top 7^ 0. Although 
the proportion of such configurations in our ensembles decreased as the temperature was increased, 
the difference Xk — Xs was still nonzero at all our temperatures. 

A study of the spectrum of the Dirac operator revealed that this breaking, at least at our highest 
temperatures (T = 190 MeV and 200 MeV), resulted from the presence of zero modes which must 
always arise whenever the underlying gauge configuration has 2top 7^ 0. Unfortunately since the 
density of the exact zero modes vanishes as V^^^^, where V is the four- volume, the observed 
breaking is unlikely to persist in the infinite-volume limit. 

On the other hand, the density of the near-zero modes is an intensive quantity. It is these modes 
that are responsible for keeping U{\)a broken even after chiral symmetry has been restored. The 
possible form of the spectral density p{X,m) that yields (yy) = but — Xs ¥^ ^ intriguing 
and currently open question. 
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